Solvable model of dissipative dynamics in the deep strong coupling regime 
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We describe the dynamics of a qubit interacting with a bosonic mode coupled to a zero- 
temperature bath in the deep strong coupling (DSC) regime. We provide an analytical solution 
for this open system dynamics in the off-resonance case of the qubit-mode interaction. Collapses 
and revivals of parity chain populations and the oscillatory behavior of the mean photon number 
are predicted. At the same time, photon number wave packets, propagating back and forth along 
parity chains, become incoherently mixed. Finally, we investigate numerically the effect of detuning 
on the validity of the analytical solution. 
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I. INTRODUCTION 

The quantum Rabi model [H |2], describing the in- 
teraction of a two-level system and a quantized single 
mode beyond the rotating- wave approximation (RWA), 
has found an experimental playground with the advent 
of novel technologies. On one hand, the development of 
solid-state semiconductors [3 and, on the other hand, 
the impressive progress on superconducting circuits [H- 
|7j, have produced the largest light-matter coupling ever 
observed. In this sense, the ultrastrong coupling (USC) 
regime happens when the coupling strength g is compara- 
ble to appreciable fractions of the oscillator frequency uj, 
0-1 ^ 9/^ ^ 1- Moreover, it is expected that these 
architectures can reach soon the deep strong coupling 
(DSC) regime, where g/oj > 1. These regimes are 
unattainable for the usual experiments in quantum op- 
tics, be in trapped ions [8j or cavity QED setups [9|.[T0]. 
In absence of dissipation, the DSC regime predicts the 
appearance of a different kind of collapses and revivals in 
the qubit level statistics, which are explained as photon 
number wave packets propagating back and forth along 
two independent parity chains [TT]. From this point of 
view, it would be of fundamental interest to study these 
key features, and possible analytical solutions, in pres- 
ence of dissipation [12]. Previous works studied related 
issues, but considering cases between the perturbative JC 
model with g/uj < 0.1 and the USC regime | T51I15) . 

In this work, we study the DSC regime of the quantum 
Rabi model in presence of mode dissipation. In Sec. II, 
we present our model and show that a zero-temperature 
Markovian bath drives the system to an incoherent mix- 
ture of two parity chains of the Hilbert space. In order 
to explore the solvability of this novel quantum open sys- 
tem, we focus on the case of a slow qubit: uiq <C {a;,^}, 
where ujq is the qubit transition frequency. This allows us 
to explain the persistence of collapses and revivals of the 
parity chain probabilities and the behavior of the mean 



photon number in the presence of dissipation. In Sec. Ill, 
we show the main results of our analytical approach. For 
the case of near-resonance, in Sec. IV, we provide a de- 
tailed numerical analysis of the asymptotic behavior of 
probability amplitudes associated to states belonging to 
both parity chains and the photon statistics. Finally, in 
Sec. V, we present our concluding remarks. 



II. THE ANALYTICAL MODEL 

We consider a general system composed of one qubit 
coupled to a single bosonic mode, as described by the 
quantum Rabi Hamiltonian 

H — huia/a + —uJodz + fig{<y + o'^)(a -I- a'), (1) 

where (Jz is a Pauli operator, a and a^ are the lowering 
and raising qubit operators, a and a^ are the annihilation 
and creation mode operators, while |e) and \g) are the 
upper and lower qubit states, correspondingly. 

It is well known that the quantum Rabi model in 
regimes where the RWA can be applied, that is the 
Jaynes-Cummings model |16j . is a solvable dynamics. 
Recently, there has been a renewed interest in analyti- 
cal and numerical approximations of the quantum Rabi 
model beyond RWA, that is, in the USC and DSC 
regimes [HI [17lf26] . Here, we consider dissipative ef- 
fects to the previous studied cases. More precisely, we 
add a reservoir acting only on the bosonic mode sub- 
system, and neglect dissipation and decoherence of the 
slow qubit. The dissipative channel will be described 
by a thermal bath at zero temperature under the Born- 
Markov approximation. This model holds for the slow 
qubit approximation and we will prove that this quan- 
tum open system still possesses analytical solutions. 

The time evolution of the above described system 
obeys the following master equation (ME) 
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where Cf[p\ is the standard Liouville superoperator 

Cf[p] — —{2apa) — a) ap — pa) a) (3) 

and K is a mode decay rate. We move into an interaction 
picture, rewriting Eq. (|l|) as H = Hq + Hi, 



Hq = hujo I -;^ + a' a 



Hi ^ hAa'^a + hg{a + a+){a + a'^), (5) 



where we have introduced the detuning parameter A = 
w — Wo, with < ujQ < uj. By means of the unitary trans- 
formation U{t) — exp{iHQt/h), we find that the structure 
of Eq. ([2| holds for the density operator pi = UpW with 
H replaced by the interaction Hamiltonian 



m = uHiW 



hAa^a + hg{e^''^°*a+a'' + a+a + H.c). 



(6) 



We consider now the qubit rotated basis, |±) — 
{\g) ± |e))/-\/2, and rewrite the total system density 
operator pi{t), obtaining the four operators Papii) = 
{a\pi{t)\l3) (a,/3 — ±), describing the bosonic mode sub- 
system. In this way, we obtain the following set of cou- 
pled differential equations 

P++ =-iA[a|at,/5++] - i.gcos(a;oi)[at + a|,/5++] 

- g su\{ujQt)[{at +a|)p_+ -/5+_(af + a\)] + Cp++ 
P — =-iA[a|at,/5 — ] + ig cos{ujot)[at +a\,p — ] 

+ g siTL{uJot)[{at +a\)p+^ - p^+{at + a\)] + Cp — 
P+- =-iA[a|at,p+_] - i5cos(wot){at -f a|,/5+_} 

- g sin{ujQt)[{at +a\)p — + p++{at + a\)] + Cp+- 

P-+ = -iA[alat,/5_+] + \g cos{uja1:){at + a|,p-+} 

+ g s\n{uJot)[{at +a\)p++ + p — {at + a\)] + Cp-+ 

(7) 

where [, ] and { , } denote the standard commutator and 
anti-commutator symbols, while dt = e~"^o*a. 

We will study how collapses and revivals behaves in 
presence of dissipative mechanisms. For this reason, in- 
stead of studying the asymptotic character of the previ- 
ous equation, we will develop an analytical method that 
works in the case of the slow qubit limit, ljq <^ {5, A}. 
We will consider evolution times under the condition 
ujot <C 1, but long enough to permite several collapses 
and revivals, wi > 1. In this case, we can approximate 
cos{u}ot) « 1 and sin(a;oi) ~ 0. 

In order to solve Eq. (l7|, we consider the charac- 
teristic function associated to the Wigner function, us- 
ing the four continuous and square-integrable functions 
X±±{a,t) = Tr[/3±±(i)£»(a)], where £>(«) = e""*""*'^ is 
the displacement operator. Here, we present the ana- 
lytical solution in the case where the qubit is prepared 



in its ground state and the mode in the vacuum state, 
/5/(0) = \g){g\ (8) |0)(0|. In this manner, we have 



X±± (a, i) = iexp{ - ^ T /3(i)a* ± P* (t)a} , (8a) 
X±T{a,t) = lF{t)cxp{ - ^ T (3{th* T P*{t)a}. 



The corresponding operators p±^±(t) take the form 



p±±{t) = -\±m){±m\ 



(8b) 

(9a) 
(9b) 



where the field coherent states amplitude /3(i) and the 
decoherence function F{t) are defined as 



/3(t)^^(c-*-l 



F{t) 



[Kt+|Qm(z*/3(t))] 



(10) 

(11) 



with the complex variable z = k/2 + iA. Once obtained 
the time evolution of the whole system density operator 
pi{t), we derive the quantities of interest such as the 
probability of the system to be in one of the states \gn) 
or I en), that is 



1 



^g/e.nV-} n'^ 



-\P{t)\- 



\m? 



l±{-lYF{t)e 



2|/3(t)l' 



(12) 

where the sign +(— ) holds for g{e). In addition, we com- 
pute the purity p[t) ~ Tr[pf (i)] of the whole system 



p{t) = 



1 



1 



:-4|/3(t)P 



(13) 



As regards the qubit and bosonic mode subsystems, we 
derive the expressions for the qubit level populations 
Pgu{t) and the photon number distribution P„(i) 



Pgu{t) = -[l±F{t)], 



Pn(t) 



-mt)\'- 



\m\ 



2n 



(14) 
(15) 



where P„ (<) is a Poissonian distribution with a mean pho- 
ton number {N{t)} = |/3(t)p. 

Given the nature of the interaction containing both 
rotating and counter-rotating terms, the Hilbert space 
of the system is split in two parity chains [Tl], associ- 
ated with the eigenvalues p = ±1 of the parity operator 
n = — (T^(— I)'' "■, that are unconnected for k — 0. For in- 
stance, the states {\g, 2N), |e, 2A^-|-1)} have parity p — 1, 
while the states {\e,2N), \g,2N + 1)} possess the oppo- 
site parity p = —1, where iV is an integer number. The 
probabilities that the system is in a state of one of the 
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FIG. 1. Effect of the parameter g/iu on probabilities Pq (t) 
(a), Pf, (i) (b) bosonic mode mean photon number {N{t)) 
(c), and purity fi{t) (d). We consider the off-resonance case 
A/o; ~ 1 and a dissipative decay rate n/uj = 0.01. 



two chains are 



Pi^\t)J-^P,Atl 



Pt\t) 



\m\ 



2n 



-Pe.oit)- 



(16a) 
(16b) 



We remark that in the unitary limit, k = 0, the decoher- 
ence function becomes F{t) ~ e^^''^'^*-" so that we obtain 
Pg o(i) = e-l/^^*)!' and Pe,o(i) = 0. Hence, starting from 
|g, 0), the evolved state vector remains in the subspace 
corresponding to parity chain p = 1. Actually, starting 
from any superposition state of the system, the time evo- 
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FIG. 2. Analytic plot of time evolution of P,i (a) and P„ 
(b) for n = (solid), 1 (dotted), 2 (dashed), (c) behaviour of 
the subsystems: mean photon number {N{t)) (solid), proba- 
bility of qubit lower state Pg{t) (dashed). For the parameters 
g/iu = 2, A/lj ~ 1 and k/uj = 0.01. 



lution of P„ (t) takes place independently in each parity 
chain [11]. Here, for k > 0, we see that Pe,oit) ¥" ^ind 
both parity chains are connected by the dissipative mech- 
anism. It is straightforward to obtain analogous results 
starting from state |e,0), just replacing F{t) by —F{t). 



III. ANALYTICAL RESULTS 

In this section, we discuss the dynamics in the limit of 
large detuning A/w ~ 1 (i.e. ujq/uj ~ 0) as described by 
our analytical model. Let us recall that this model is only 
valid for dimensionless time t = tot <^ 1/(1 — A/w). The 
focus of this discussion will be on the peculiar dynamics 
of the DSC regime tH], such as collapses, revivals, and 
parity chains and how they are modified by the presence 
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FIG. 3. Analytic plot of probabilities Pi*'(i) for n < 20 
sampled at times n — nl (with Z = 0, . . . , 12), with the same 
parameters as in Fig. [2] 



of dissipation, k. 

Figure [T] shows how the DSC dynamics emerges under 
the condition A/w ~ 1, as one increases g/uj- The first 
signatures are the appearance of collapses and revivals in 
the probability Pg {t), and a periodic dynamics of the 
average photon number {N(t)), as it was already seen 
in Ref. [TT]. The novel feature consists on the transfer 
of probability from the positive to the negative parity 
chain, evidenced in Fig. [i]d. This effect is induced by the 
presence of the small dissipation, n = O.Olw. Another 
consequence is the fact that the purity /i(t) decreases by 
increasing g/uj and that, in the DSC regime, the system 
reaches quite rapidly a maximally mixed state. 

In Fig. [2] we show with greater detail the time evo- 
lution of the two parity chains, setting g = 2a;. Proba- 
bilities P„ and Pn exhibit revivals at r = 27r x Z"*", 
and there is a progressive transfer of it from the positive 
to the negative parity chain. The consequence is that, 
at the end, the revivals of both chains exhibit the same 
structure and height, and that the qubit ends up with 
the same probability in the state \g) or |e) (cf. Fig. [2];). 

In order to better illustrate the role of dissipation in 
connecting the two parity chains, we plot in Fig. [3] differ- 
ent time shots of the statistics P„ (t) for large enough 
values of n. We notice that, starting from \g,0), the 
two chains are connected via the dissipative channel re- 
sulting in an incoherent mixture of them, breaking their 
independent dynamics stemming from the unitary case. 

In our analytical framework, we can compute the 
steady state of the system. 



PS 



\(3s){(3sm+){+\ + \-Ps){~Psm-){-\ , (17) 



where Ps = —^g/z is the steady amplitude of the field 
coherent state. The structure of the steady state is note- 
worthy, since it is a statistical mixture of two parts, one 
that associates j/Ss) to the qubit state |-|-) and the other 
connects | — Ps) to |— ). 

This means that the system relaxes upon the two 
states |±, i/Ss), which have the same energy mean value 
Ap, = h/\\l3s\'' + hg{ps + P*s)- In Fig. g we show 
how the first elements of the two parity chains P„ (t) 
(n = 0,1,2) and the mean photon number {N(t)) ap- 
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FIG. 4. Analytic plots of the approach to the steady state of 
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(i) (a) and PA ' (t) (b) for n = (dashed), 1 (dotted), 
2 (solid), (c) Mean photon number {N{t)}, with g/cj = 2, 
A/o; ~ 1, and k/cj = 0.2. 
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FIG. 5. Comparison between the mean photon number 
{N{t)) from the analytical solution, in the limit A/a; ~ 1 
(sohd), and through MCWF method for A/o; = 0.8 (dashed), 
in the DSC regime g/u — 2 and for k/cu — 0.01. 
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FIG. 6. Probabilities P^'^^t) (solid) and P^^''' (t) (dotted) as 
a function of A/a; in the DSC regime, g/cj = 2 and k/uj = 
0.01. 



proach the steady state for relatively large values of de- 
cay rate k/lo. We remark that Pq ' (t) vanish at long 
times, while the other probabilities reach a non-zero con- 
stant value, so that the steady state Poissonian photon 
statistics P„(i) = Pn (t) + PrT (t) has a mean value 
|/3s|2 = 4(7V(k^ + 4A2), with A/w ~ 1. 
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FIG. 7. Effect of detuning in the DSC regime g/ui = 2 
and k/lu = 0.01. Plots of {N{t)): analytical (solid-gray) and 
numerical solutions for A/cj = 0.75 (dashed), A/uj = 0.5 
(dotted), A/uj = 0.25 (dash-dotted), A/co = (solid-black). 



IV. NUMERICAL RESULTS FOR < A/lj < 1 

The role of this section is to investigate numerically all 
regimes of the detuning A using the Monte Carlo wave- 
function method [27]. As a first example, we shall dis- 
cuss a set of simulations for small values of the detuning 
A/w = 0.8, for which the analytical solution may still be 
valid at short times t ^ 5. These are values that may 
be achieved for experimental circuit QED setups [51 [5] , 
selecting a resonator of a; ~ 6 GHz and ojq — 1.2 GHz. 
The results of the simulation are shown in Fig. [5] where 
we plot the Monte Carlo simulation together with the 
analytical ansatz. It is evident that both solutions are 
very close, with an error < 4% for r < 9, justifying the 
approximations of Sec. II, at least for 0.8 < A/w < 1. 

We investigate now the effect of moving towards the 
full resonance A/uj ~ 0. The behavior of parity chain 

probabilities P„ (t) is heavily affected by the reduction 
of the detuning values A/w. As an example, in Fig. |6J 

we show that the probabilities Pg (t) are progressively 
spoiled, losing the symmetry of the ideal case A/ui ~ 1 

(see Fig. [Sp,) . We outhne that the probabilities P„ ' (t) 
with n > present irregular oscillations that do not van- 
ish for long times. Despite the distortions of P„ '{t), we 
show in Fig. It] that the mean photon number {N{t)) ex- 
hibits regular oscillations. The only effect of decreasing 
the rate A/w, towards a resonant condition, is to reduce 
the amplitude of oscillations and to induce a small shift 
in the peak times. In order to quantify this effect, we 
estimate the error in each curve at r = 8.5 with respect 
to the analytical solution. The percentage of relative dif- 
ferences are 7% (A/w = 0.75), 17% (0.5), 29% (0.25) and 
38% (0). These features can be explained as an internal 
dephasing mechanism due to the qubit free energy [llj . 
Finally, let us consider the simultaneous effect of differ- 
ent decay rates and detunings. Figure [8] shows the results 
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FIG. 8. Mean photon number {N(t)} for difTerent values of 
parameter k/lj: 0.3 (solid), 0.5(dashed), l(dotted) and for the 
coupling parameter g/co = 2. (a) Numerical results composed 
by four curves for each k/uj value, corresponding, from top to 
bottom, to A/lo = 0.75, 0.5, 0.25, 0, (b) Analytical solution in 
the far off resonant limit A/lj = 1. 



of the numerical simulation for the mean photon num- 
ber. We notice that the effect of a decreasing detuning 
becomes almost negligible in the presence of dissipation. 
Moreover, the asymptotic value of {N{t)) agree very well 
with the analytical predictions of Sec. II. 



V. CONCLUSIONS 

We have studied the DSC regime of the quantum 
Rabi model with a zero-temperature bath interacting 
with the single quantized mode. We have found analyt- 
ical solutions for the off-resonant case, describing the 
dissipative dynamics that induces incoherent mixtures 
of both parity chains. Furthermore, we have presented 
a numerical analysis for the near-resonant case where 
analytical solutions are not available. Finally, we have 
estimated the limits of our analytical solutions as A/w 
decreases, among other features. 
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